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Motivated by understanding the emergence of thermodynamic restoring forces and oscillations, we 
develop a quantum-mechanical model of a bath of spins coupled to the elasticity of a material. We 
show our model reproduces the behavior of a variety of entropic springs while enabling investigation 
of non-equilibrium resonator states in the quantum domain. We find our model emerges naturally 
in disordered elastic media such as glasses, and is an additional, expected effect in systems with 
anomalous specific heat and 1/f noise at low temperatures due to two-level systems that fluctuate. 


I. INTRODUCTION 


The resonant, linear response of a variety of physical 
systems naturally leads to using harmonic oscillator ap¬ 
proximations to describe their behavior. Recent break¬ 
throughs in the fabrication and characterization of large 
mechanical oscillators have led to tremendous progress in 
exploring the quantum behavior of macroscopic systemsi, 
including reports of ground state cooling^Ti^ and demon¬ 
strations of squeezed states beyond the standard quan¬ 
tum limit^. However, not all that resonates is a quantum 
harmonic oscillator. A variety of materials exhibit elas¬ 
tic, Hooke’s law-type behavior at high temperatures due 
to changes in their microscopic configurations, i.e., their 
entropy2r— . A typical example of such “entropic springs” 
is given by a rubber band, where forces applied cause dis¬ 
ordered strands to straighten, reducing their entropy and 
resulting in a restoring force near thermodynamic equi¬ 
librium. 


Here we show how a resonator whose stiffness de¬ 
pends upon the configuration of rapidly thermalizing 
two-level systems (TLSs) is a natural example of an en¬ 
tropic spring. This type of system arises in a variety of 
atomic and condensed matter systemsi^^— . Specifically, 
fast thermalization of the TLSs in such systems leads 
to the first quantum-mechanical model of an entropic 
spring. The corresponding quantum model is surprisingly 
simple and admits solutions in a variety of physically rel¬ 
evant scenarios. We consider a collection of TLSs coupled 
to the elasticity of a medium, and find this naturally sat¬ 
isfies the fast thermalization mechanism and results in 
an entropic spring restoring force. Furthermore, entropic 
oscillations persist in the thermodynamic limit at any 
finite temperature, and we show numerically that in typ¬ 
ical cryogenic environments mechanical resonators used 
in optomechanics have a non-vanishing entropic spring 
contribution. Nonetheless, quantum coherence persists 
in the thermodynamic limit. 


II. QUADRATIC-COUPLED SPIN BATH 
MODEL 


As a simplified model, consider a ID harmonic oscilla¬ 
tor with a mass M and a bare frequency ujq interacting 
quadratically with a collection of independent TLSs. The 
Hamiltonian of the combined system is Us = LJq+E with 


i/n = 


Mloqx'^ 


2M 


N 


( 1 ) 


where the TLSs are described as pseudo spin-1/2 systems 
using Pauli matrices cr/, N is the number of indepen¬ 
dent TLSs, and B is the energy difference of the TLSs. 
Motivated by our analysis of the TLS defects below, we 
consider a coupling term between the oscillator and TLSs 
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We denote this the quadratic-coupled spin bath model 
(QCSB), which is a also simplified form of strong coupling 
Tavis-Cummings model in dispersive regioni^. It can be 
pictured as a system whose effective spring constant is de¬ 
termined by the polarization (a) = ((cr* + l)/2)/N. 

We are primarily interested in the thermodynamic limit 
A —^ oo. Here we consider homogeneous case for simplic¬ 
ity; distributions of parameters will be considered later. 

The thermodynamic properties of the system are en¬ 
coded in its partition function. Specifically, the quantum 
partition function of the QCSB model can be evaluated 
as Z = Tre“^^, /3 = l/{kBT) is the inverse temperature. 
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= ^TLS^OSC- (3) 


We drop terms of 0{^) in the thermodynamic limit us¬ 
ing the method of steepest descent, is the parti¬ 

tion function of N independent TLSs, and Zqsc can be 
interpreted as the quantum partition function of an os- 
ci llator with a temperature-dependent frequency WeS = 
y/ojQ + (o’)^, and polarization {a)g =1/(1 + e^^). 
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We proceed to see whether the system responds to an 
external perturbation as an entropic spring. We add 
an external field Hf = —fx corresponding to a con¬ 
stant force /. The partition function becomes Z ss 
Zo exp(/3/^/2Ma;gfj) for small /. 

The Helmholtz free energy under the perturbation is 
A = U - ST = -\n{Z)ll3 = Aq - p/{2Muj^g), with 
unperturbed Helmholtz free energy Aq = — In(Zo)//?. 
Crucially, a portion of the change in A arises from the 
change in entropy of the TLSs, S = —^ = Sq — 




Qn — _£do 
<50 — Qx ■ 

This entropy change arises microscopically by assum¬ 
ing that the spins thermalize much faster than the mo¬ 
tion of the oscillator. Specifically, consider x now as a 
constant. Under fast thermalization, each spin goes to a 
Gibbs state according to the instantaneous energy differ¬ 
ence as 


J _ 1 o^^ 

^TLS j 0 l) ' 


( 4 ) 


Here 5 = MuP/2N. The total entropy of the spins is thus 
a function of x expressed as 


S'tls ^ Sb — ks 


2(l + e/5S)2 ' 


( 5 ) 


where 5 'b is the total entropy of spins with energy dif¬ 
ference B. Adding a force / to the system will dis¬ 
place the oscillator equilibrium position by an amount 
Ax = f/{Mpg) and decrease its entropy by AS = 
u^pMuZBpp^ArV - kK The ore 

viously found entropy change is recovered. 

From the first law of thermodynamics and eq.®, there 
will be a corresponding static entropic force 


We now derive the entropic spring master equation un¬ 
der the conditions of fast spin thermalization. The below 
derivation follows a standard quantum optics approach!^. 
Taking the total Hamiltonian H = Hs + Hr + Hsn, we 
have the Hamiltonian of the oscillator-spins system in the 
QCSB model: 
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the Hamiltonian of the thermal reservoir consisting of 
independent bosonic modes ^ ^ with frequencies ui^ ^ ^ 
for each spin j: 

= E Ai^A ’"feAi; (9) 

fe,A., ” ” 

and the interaction between the spins and the reservoir: 
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( 10 ) 




We remark that the spin-reservoir interaction is under 
the dipole approximation and the rotating wave approxi¬ 
mation, appropriate for B A , and the reservoir only 
interacts with the spins. 

To focus on the system-reservoir interaction, we now 
change to the interaction picture. An operator in the 
interaction picture 0{t) is defined through the unitary 
transformation 


Fs = T 


SStls 

dx 


MPpBe^^ 
{l + ePBf 


X = Mujgx. 


( 6 ) 


Thus, a portion of Hooke’s law arises entirely from re¬ 
configuration of the TLSs. We note that this isothermal 
process is reversible. We define an entropic spring pa¬ 
rameter Rs to quantify the entropic contribution, 


Rs = 




ujpPBe^^/{I A 

^eff 


( 7 ) 


the ratio between the frequency associated with entropic 
Hooke’s law (eq. IH]) and the effective frequency. 


d{t) = ^ 


( 11 ) 


and the time derivative of the interaction operator be¬ 
comes 


6{t) 


HsR{t), d{t) 


( 12 ) 


Defining = ujq S ^ , the oscillator opera¬ 

tors in the interaction picture can be solved as follows: 

P 

xit) = xcos{uJY,t) -\ -sinlwst), 

Mwt. 

p{t) =pcos{u)Y.t) — Mujsxsin{u}st). (13) 


III. ENTROPIC SPRING MASTER EQUATION 

The proposed QCSB model behaves like a spring in 
static equilibrium. However, without thermalization of 
the spins, the model will not undergo entropic oscilla¬ 
tions. Thus, a full quantum treatment necessarily in¬ 
volves the thermalization of the spins. 


Since we are working in the regime where spins will 
quickly thermalize to their equilibrium value at temper¬ 
ature T, ~ TTTra- + P = I/UbT, we can 

approximate x(t) ~ x{t) = xcos{uJa-t) + sipLUg-t) 
using a spin-independent frequency w? = ojq P/(\ + 
when evaluating spin evolution (while the coupling 
strength between x and each spin is already of the order 





















3 


of 5). We now examine the spins and the bath in the 
interaction picture: 


Starting from the exact time evolution of the S ® R 
density matrix in interaction picture 


a+{t) = 

k,X,j k,X,j ’ 

k,X,j 

^ ■ (1"^) 
k,\,j 

Let x(t) be the density operator for S ® R, then 
p{t) = tr/j [x(t)] is the reduced density matrix describ¬ 
ing the system only. Assume the interaction is turned 
on at t=0 and the total Hilberspace starts with an un¬ 
correlated (product) state x(0) = x(0) = p(0) ® Rq. 
The intial state of the reservoir is taken to be a thermal 
equilibrium state at temperature T, with density matrix 

fcAd_I 


X{t) = -i Hsr,x{0) 


- / dt' 
Jo 


HsR{t), HsR{t'),x{t') 

(15) 


Assuming weak coupling between the reservoir and the 
spins, we first introduce {i)Born approximation, so that 
x(t) ~ p{t) ® Rq. The time evolution of p{t) = Rr [x(i)] 
can then be evaluated through the Born approximation 
as 


p(i) = -y dt'trR^^HsRit), HsR{t),p{t')Ro 


(16) 


Plugged in the interaction Hamiltonian for our model, 
the time evolution equation eq. (1161) can be explictly writ¬ 
ten down as 


^(t) = -J 2 [i?or](t)f,(t') 


h.c. 




i?orUi)f,(t') 
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tVR 


R^T,{t)T]{t') 


+ h.c. 
h.c. 




Assuming short time correlation between the reser¬ 
voir and the spins, we now make the second approx¬ 
imation, (ii) Markov approximation, so that operators 
change slowly with time within the integral. We can then 
replace p{t') by p{t), qr{e-^s j^xRt,)dt,^ 

T{e~'‘^ Jo a;^(ti)dti| ^y ^-iSx{t)t' formula, and 

extend the lower bound of the time integral from 0 to 

— CXD. 


After applying the above approximations, terms with 
di. on the same side of p{t) resemble that of damped 
two-level atoms with the thermalization rate associated 
with a frequency B + 5x^{t) depending on x. Terms with 
(7^, on opposite sides of p{t) can be expanded through 
a complete basis of x{t) for further calculation. Take the 
a^^p(t)ai_ term for example. 




' —OO 
rt 


f "7 


dt' / dxidxa |a:i(t)) (a:i(t)| cr^p(t)ai |a:3(t)) (a:3(t)| t i?or](t)l 7 (^') 


' —OO 
^OO 


dr / dxidxs |xi(t)) (xi(t)l 4p(t)ai |x3(t)) (x3(t)| [i?ori(t)r;(t - r)! . (18) 


In the last equality we perform a change of variable using t = t — t'. 
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Plugging in the expressions of reservoir correlation 
functions and evaluate the integral through standard 
master equation derivation approach, the time evolution 
equation becomes 

P{t) = {<cri 7 $(t)(«$(t) + l)p (0 

3 

+ Pit)ins(t) + l)7$(t)0'+0'- 

+ cr^-(J+lx{t)nx(t)P{t) + P{t)nx(t)lx(t)(^-(^+ 


(a) (b) 



Here we have 7^(4) = 2 ttJ{B + Sx{t)), the spin ther- 
malization rate associated with energy B + Sx(t), and 

nx(t) = ile-i 3 (B+ 6 x(t )) ; the thermal occupation number 
associated with energy B + Sx{t). x{t) is the interaction 
position operator and J{u!) is the spectral density of the 
reservoir. We have taken Lamb-shift terms to be zero by 
appropriate redefinition of B. 

Now we transform the time evolution equation back to 
the Schrodinger picture using 

m = -i Ws, Pit)] + (20) 

The complicated time-dependence in ea. dT^ disap¬ 
pears neatly when transformed back to the Schrodinger 
representation. When the thermalization rate is faster 
than the mechanical oscillation, we find that the dynam¬ 
ical evolution of the entropic spring model can be de¬ 
scribed by an entropic spring master equation 

Pit) = -i \Hs,pit)] - ^ {o-:(_cri7a;2(h,r2 + l)p(t) 

3 

+ pit)inx'^ + l)7x2cr^cr(_ - a^_-fx2inx2 + l)p(t)cr+ 

- CjLpit)inx2 + l)7a;2Cr+ -b cricr+ 7 a; 2 ha; 2 p(t) 

-b pit)nx2^x'2cr^-cr+ -cr+ 7 x 2 ha; 2 p(t)(Ti - a^^pit)nx2-/x2'ji_ 

( 21 ) 

Here we have the thermalization rate 72,2 = 27r J(i3 -b 
Sx^), and the thermal occupation number 71 , 3,2 = 

^-D(B + 5x'2) 

1 — q — ^(B + Sx2) ■ 

Note that the above master equation decouples to a 
Lindblad master equation for damped TLSs and the uni¬ 
tary evolution for a harmonic oscillator when taking the 
coupling w to zero. In the massive limit, in which the 
oscillator position can be treated as a static constant, 
the spins will thermalize according to the probability 
= e-P(B+5x ^vvhich is consistent with the detailed 
balance as suggested by a static equilibrium. 


FIG. 1. Numerical simulation based on the entropic spring 
master equation with parameters loq— 0 , cj = Vl -b 7 = 
10, /3= 0.01, B — 100, and several N as identified in the 
legend, r is the oscillation period r = 2'K/iJeS- [h = M = 

UJeS ~ 1 ) 


Observables associated solely with the oscillator, i.e., 
operators which are functions of x and p, will evolve ac¬ 
cording to the master equation eiD as 

(oix,p, t)^ = Tr {Opit)} = -i ([O, iLs]), (22) 

using the cyclic property of trace and the fact that x, 
and p commute with spins. The overall contribution from 
the dissipation part vanishes. That is, (x) and (p) evolve 
only through the Hamiltonian. In the thermodynamic 
limit, we recover our static result: the oscillator evolves 
as a simple harmonic oscillator with an angular frequency 

Weff = y'ujQ +U}'^ {^)b- 


IV. MESOSCOPIC EFFECTS 



FIG. 2. (a) Inverse of position oscillation damping rate versus 
N. (b) Inverse of dephasing rate versus N. (c) Dephasing rate 
versus different entropic spring contribution. We use the same 
set of parameters as Fig. [T] for (a)(b), and use parameters 
iV = 10, w = {0.7,0.71,0.72, ,1}, ujo = 7 = 10, 

/3= 0.01, B = 100 for (c). r is the oscillation period r = 
27 r/tUeff . {h = M = OJeff = 1) 
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We now consider mesoscopic corrections. The entropic 
spring master equation (eqUU cannot be solved analyt¬ 
ically due to nontrivial correlation between x and cr| so 
we turn to numerical simulations. To overcome exponen¬ 
tially large density matrix representations, we use sym¬ 
metric properties of the system Hamiltonian to reduce 
computational degrees of freedom. We define a set of 
elements {C^}, j=0, 1, ..., N, symmetric in spins 

N N 

Co = I,C,=Y,al,C2=Y,Y.a>i ,••• , 

2=1 2=1 j >2 

Cn = . (23) 

All completely symmetric density matrices of N spins can 
be expanded as Ps = then construct 

the transformation laws for '^jO'iCk, '^j<J+CkCfl_ and 

'^ja'^_Ckcr\. and reduce the spin dimension from 2 ^ to 
N-i-1 with the above algebra. 

Our numerical simulations (FiglT]) follow the quantum 
trajectory of an initial state chosen as the Gibbs state 
of N spins with energy difference B and a coherent state 
I a) with a = \ for the oscillator. For our simulation, 
we use units such that h = M = uJeS = 1. We take 
J{uj) ~ constant near B for simplicity, and set 7 large 
compared with ujes to achieve fast thermalization and 
focus on cases with wq = 0 corresponding to oscillations 
arising entirely from reconfiguration of the coupled TLSs; 
cases with wq 7 ^ 0 behave similarly. 

Figure [T] shows a simulation of the entropic spring evo¬ 
lution for several N. In Fig. [IKa)(b) we see damped os¬ 
cillations in {x{t)) and (x^(t)) with damping rate sup¬ 
pressed as N increases. Note that in contrast with 
eq. ( 1 ^^ . we see damping in oscillations due to uncer¬ 
tainties in the oscillation frequency caused by finite N. 
The entropic nature of the system is described by the 
TLSs von Neumann entropy Stls = —Tr[pTLS ln(/OTLs)], 
where ptls = Ti'oscP is the reduced density matrix de¬ 
scribing the TLSs. Fig (TJc) shows the evolution of 
the TLSs entropy minus its initial thermal equilibrium 
value, AS'tls = •S'tls ~ 5'thermai, which oscillates with 
— (a;^(<)) through an oscillating amplitude independent 
of N. These results agree with our requirement that 
the spins are constantly rethermalizing according to the 
value of x^, and S'tls is changing with x^ with a N- 
independent coefficient as suggested by eq.®. Finally 
we monitor the dephasing by plotting the impurity of 
the oscillator, l-Tr[pQgj,(t)], in Fig. [IKd), where Posdt) = 
TrTLs[p(i)] is the oscillator reduced density matrix. The 
dephasing reduces as N increases. 

The qualitative behavior of the damping of {x) and 
dephasing of l-Tr[yOQgj,(t)] is followed in Fig. [2l By fitting 
the evolution with an exponential decay envelope, the 
damping rate of {x(t)) is found to be proportional to 1/N 
(Fig. [IKa)). The dephasing is also suppressed by 1/N as 
shown in Fig. [Hb) and is proportional to the entropic 
spring parameter Rs as shown in Fig. ®c). 

The numerical simulation shows the intriguing result 


that, in thermodynamic limit, the dephasing rate of 
the oscillator vanishes while amplitude of the oscilla¬ 
tion of the entropy 5'tls stays constant. This phe¬ 
nomenon can be explained qualitatively by interpret¬ 
ing the thermal spins as a continuous weak measure¬ 
ment of x'^ through the coupling. The flipping rate of 
TLSs is Fflip = 7 x 2 (h:j ,2 -h l)P(t) + “ ^(t)) ~ 

272 , 2 e“^(^+‘^“^V(l assuming 73,2 is a con¬ 

stant and under fast thermalization approximation. The 
measurement rate is proportional to the flipping rate con¬ 
dition upon the change in x^, which is of the order of 
a;^ 9 a; 2 rflip (X ^(1)) causes the N-independent 

entropy oscillation amplitude. The measurement back- 
action depends on x"*, therefore the dephasing rate is 
proportional to ~ C>(l/fV), which vanishes in 

thermodynamic limit. 

In principle, the high temperature limit of this model 
should be treatable in a semi-classical approximation 
where the quantum nature of the spin bath is largely ig¬ 
nored. However, such treatments would necessarily miss 
the mesoscopic dephasing effects due to the effective back 
action of the spin bath on the resonator shown in Fig. [1] 
and Fig. [2] 


V. MICROSCOPIC ORIGINS: TWO LEVEL 
SYSTEMS IN AMORPHOUS SOLIDS 


We now consider physical systems described by our 
model. Our QCSB entropic spring Hamiltonian follows 
naturally from the interaction between a single phonon 
mode and a collection of defects in amorphous solids. De¬ 
fects in such materials can be microscopically modeled as 
atoms which can tunnel between two local ground states 
of asymmetric double-well potential a^^d^ . The effective 
Hamiltonian of a TLS defect is 


n ^ ^ 

xItLS = - —ax 


1 

2 



(24) 


where A is the asymmetry of the double-well, and A is 
the tunneling between the wells. This model is used to 
explain anomalous low temperature specific heat behav¬ 
ior in glassy solids^^, and yields a I /f noise spectrum in 
macroscopic properties such as electrical resistances^. 

An external strain field adds a perturbation to TLS 
Hamiltonian of the formSS 

Hp = - '^lijSij fg (25) 

ij ^ ^ 


where s = lij are the coupling constants, 

and Sij is the strain field tensor. The displacement field 
caused by a single longitudinal wave can be expressed 
as quantized phonon creation and annihilation operators 

^ = \J 2 J^Vp (®fc + afc) cos{kz)z, where ujk = VLk, 
k is the wave vector and vl is the sound velocity for 
longitudinal waves traveling in the solid. The strain field 
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displacement u(x) 



FIG. 3. A single longitudinal phonon mode interacts with a 
collection of two level system defects in a solid. 


T(K) 

1.4 

1.2 
1.0 
0.8 
0.6 
0.4 
0.2 

20 40 60 80 100 

cOp (MHz) 



Rs 

0.0009 

0.0008 

0.0007 

0.0006 

- 0.0005 

- 0.0004 

- 0.0003 
0.0002 
0.0001 


tensor is followed by Sij = (djUi + diUj)l2. Assuming 
the interaction is isotropic, the coupling term is of the 

form Hp = go'z^ 21 ^^ + Ofc) with coupling strength 

g. For transverse waves the above analysis holds with 
replaced by -jt, and vl replaced by vt- 

Changing to the diagonal basis of Htls^ the coupling 
between a single phonon mode with frequency ojp and 
intrinsic TLS defects inside is of a general form 



Mujpx"^ 

2 



(26) 


Here we have x = ^ 2 Muip — a), 

tj = is the excitation energy for the TLS, 

gx,j = -ffAj/cj, and = g^jl^j- 

We are working in the physical region where the TLS 
energy differences e is the dominant energy scale in the 
system, Wp, 5 ^, 5z ^ £• We first rotate the spin 
about (T^ with a small angle 9^ = xgxj/{xgz,j + ej/2) to 
diagonalize ejai/2 + x(gxj<jl + g^jai), which also boosts 
p to p — (ytgx,j/i ^3 + The generated a-j. term in p 

3 

is negligible in a rotating-wave approximation since €j is 
the largest energy scale in the system. 

We then perform a polaron-type transformation 
U = to remove terms linear in x, in order 

to better approximate the dressed vacuum, with ^ = 

-~ J^^z assuming the fluctuating 

TLSs have a minor effect on existing mechanical oscilla¬ 
tion. 

Finally, we have 


M(a;^ - ^V2) _ Ma;| aj + 1 . 

2M ^ 2 ^22 


+E 


j 



E 


gz,igz,j 

2MujI 


(27) 


Here a;| = Yhj , M = Vp. Note that the 

spring constant contribution from each TLS is oc 1/M oc 


FIG. 4. Entropic spring parameter (Eq. 1281) as a function 
of resonator frequency cUp and bulk temperature T for pa¬ 
rameters in the text corresponding to a BAW resonator in 
Si02 similar to Ref.—. At intermediate temperatures around 
a tenth of a percent of oscillation is entropic in nature. 


1/N. Our entropic spring Hamiltonian now arises natu¬ 
rally from this physical system, up to an extra Ising term. 
The strength of the Ising-like term is negligible compared 
to Cj, and is thus omitted in the thermodynamic analysis. 

Based on the above, we can estimate the entropic 
spring parameter. Unlike the uniform case, as our ini¬ 
tial QCSB model, we have to sum over the distribution 
of TLSs physical parameters. The entropic spring pa¬ 
rameter Rg is 


Rs = 


— 


d (l+e'3'5) 


\ < - T., 


^9l 


J A#ltanh^ 


(28) 


For a standard tunneling TLS model that repro¬ 
duces experimental 1/f noise, the probability distri¬ 
bution function for a TLS with energy difference e 
and relaxation rate F can be expressed as P(e,r) = 
P/(2rv'l-r/F^ax(e))^. Here we have a constant 
density of states P, and the TLS relaxation rate P 

as r(e,A) = (|+2|)5ig,coth(/3e/2)^^. The 
maximum relaxation rate for an energy difference e is 
rmax(e) = r(e, e) since A < e. 

The summation over TLSs is equivalent to integration 
over probability distribution, V f de f drP(e,r). 

At temperature T, the contributing TLSs are those of 
e < ksT, since the entropic spring model works in the 
physical region that the TLSs fluctuate, which sets Cmax 
to ksT. To be consistent with previous analysis, we con¬ 
sider F > Wp, which sets the lower limit of F to ujp and 
Cmin such that rmax(Cmin) — ^p- 

By assuming small displacements, we approximate the 
sinusoidal displacement field to a triangle wave, and the 
variation of cos{kx) in x is now a constant slope A up 
to a ± sign. For a phonon polarization denoted by a, 
using the fact that A^/e^ = T/Tmaxp), Pxj is identified 
as glj = jlkxjpT/Air'^vlVpTjnaxp)- The entropic spring 
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parameter can be evaluated through the integration over 
the above identified parameters. 


There has been recent success in trapping longitudi¬ 
nal phonons in quartz Bulk Acoustic Wave (BAW) res¬ 
onators with high quality factors^^. For the case of Si02, 
we plot the lower bound of Rs as a function of longitudi¬ 
nal mechanical frequency ujp and temperature T (FigH]), 
using the parameters given in reference^. The maximum 
value of Rs is about 0.1%, and there is no entropic con¬ 
tribution for rmax(A:BT) < Wp, which happens when we 
have 2.86 x < Wp, and can be identified by 

the region below the red line in Fig. SI 


VI. OUTLOOK 

Our quantum entropic spring model predicts non¬ 
trivial contributions to the spring constant in disordered 
systems. However, direct observation of dephasing in¬ 
duced by the effect remains challenging. Future efforts 
in synthetic versions of this model, such as in the Tavis- 
Cummings interaction of an ensemble of atoms in an op¬ 
tical cavity, may allow experimental observation of the 
effects if an appropriate fast thermalization mechanism 
of the atoms can be provided. 
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